00 

o 



Response Function of an Irregular Oscillator 

Hirokazu Aiba^ , and Toru Suzuki^ 

^Koka Women's College, 38 Kadono-cho Nishikyogoku, Ukyo-ku, 615-0882 Kyoto, Japan 

^Department of Physics, Tokyo Metropolitan University, 192-0397 Hachioji, Japan 

(February 5, 2008) 

Properties of the response functions for a two-dimensional quartic oscillator are studied based 
on the diagonalization of the Hamiltonian in a large model space. In particular, response functions 
corresponding to a given momentum transfer are studied for different values of the coupling param- 
eter in the Hamiltonian. The latter controls regular or chaotic nature of the spectra and eigenstates 
of the system. Fluctuation properties of the energy-strength correlation of the response are inves- 
tigated. Even when the statistical properties of the system indicate an almost completely chaotic 
character, there remains a typical structure in the response function similar to that in the regular 
system. The nature of this structure is studied in some detail. 
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Quantum mechanical manifestations of dynamical properties of a system which classically possesses a chaotic 

CN ' character have been intensively studied KM . Level statistics which has a long history in nuclear physics as described 

^ ! ^y ^^^ random-matrix theory ^j^ is now a favorable playground in the discussion of a transition from a regular 

(integrable) to a chaotic character of a quantum system. Together with numerical studies on model systems, analytic 

investigation has been made based on semiclassical trace formula ||] . Wave functions of a system which is classically 

p^ , chaotic have also been investigated: Statistical theory predicts that the amplitude distributions show the Porter- 

^—{ • Thomas distribution Q, which was then numerically demonstrated to hold for model chaotic systems. Contrary to 

0^ the naive expectation, however, the profile of the wave function for a chaotic system is not entirely structureless: For 

0^ ■ instance, the Husimi representation of a wave function in the chaotic regime frequently suffers from a scar of classical 

r^ [ periodic orbits 0]. Although a considerable progress has since been made, it is still an important issue to clarify the 

J>^i characteristic of eigenstates and its matrix elements for systems which are classically irregular or chaotic. 

''O ' It is the purpose of the present paper to study another aspect of the wave functions for a system which shows 

J^^ , a transition from an integrable to a chaotic character: We study response functions of the system, i.e. transition 

jrt ' matrix elements as a function of energy. Statistical properties of the distribution of transition matrix elements have 



^ 



f~| , been studied |pHlO|, and it was shown, in particular, that the distribution becomes a Porter-Thomas type for chaotic 
O ' systems. The approach proposed in Ref. ||l^ has since been developed to elucidate the role of periodic orbits and was 
^ ■ extended to various systems including the response of mesoscopic systems to realistic probes pl|-p^ . These studies 
are based on the semiclassical framework, and focus mainly on the responses to long- wavelength probes. Semiclassical 
studies of response functions have also been done in a different framework p5[ | , which concentrate on their smooth 
behavior but not much on fluctuations. 

In this paper, we study a model system using a large space diagonalization and calculating response functions for 
the operators which probe the system with a variable wavelength, or momentum transfer. In particular, we put an 
emphasis on which aspect of the response function reflects the regular or the chaotic character of the system. We 
also would like to study a structure in the response function for a system in a chaotic regime which however is not 
expected to occur in the statistical random-matrix theory. Although we are here concerned with the properties of 
response functions which show up in a model system, they will also be interesting in realistic applications, as this 
response is similar to the excitation cross section for, e.g., electron scattering in the plane- wave Born approximation. 
Thus it is hoped that the present study may provide insight into the understanding of collective states such as nuclear 
giant resonances embedded in complicated many-particle many- hole states as studied in nuclear reactions. 

The main content of the paper is as follows: In the next section we summarize classical and quantum mechanical 



properties of the model Hamiltonian and flx values of the relevant parameters. In Sec. Ill we study response functions 
first for a long wavelength probe, and then for the probe characterized by a given momentum transfer. We study the 
fluctuation properties of the response functions, concentrating especially on the similarity or the difference for the 
regular and the chaotic systems. The accuracy of the calculation has been checked against sum rules. Final section 
is devoted to a summary. 



II. BASIC INGREDIENTS OF THE MODEL 

In order to study the response functions for a system which is capable of showing regular as well as chaotic properties, 
we adopt the following Hamiltonian as a model, 

H =\{pI+pI) + \{x^ +y^)~kx^y\ (1) 

This model Hamiltonian has been adopted by a number of authors for the studies of level statistics or wave functions 
|l7|-pi|. It was also employed as a model for a background system in the studies of the fluctuation properties of 
strength functions. [p3| Let us first briefly summarize classical properties of the model, which have been studied in 
detail by Meyer iQ . The Hamiltonian (p possesses a dynamical scaling property in the sense that the classical phase 
space structure at one energy is mapped into another by a simple scaling of the coordinates and momenta. It has a 
high symmetry called Ci^: The Hamiltonian is invariant with respect to a reflection about a;-axis, y-axis, and also 
about the line x = y. Furthermore, by rotating 45° in the x — y plane, the Hamiltonian is mapped into the one with 
a coupling constant k' = —(3 + k)/{l — k). As the system becomes unbounded from below for fc > 1, we have only to 
consider the range [—1,1] for the coupling constant k. Meyer ||l8| showed that for large k values (> 0.4) the classical 
phase space structure is almost completely chaotic, while for small k the system becomes regular. In the following 
calculations we adopt two typical values of the parameter: k =0.2 and 0.6. They correspond to quasiintegrable and 
fully chaotic systems, respectively. For instance, for k — 0.6 a single trajectory fills up almost 90% of the available 
phase space, while for k = 0.2 the fraction of the phase space covered by irregular orbits is only 25% in typical cases 

®- ... . . ^ 

To study quantum mechanical properties of the eigenstates of the Hamiltonian (|lj), we follow the procedure es- 
sentially of Zimmerman et al. p9[ : The Hamiltonian is diagonalized within a truncated model space spanned by a 
set of suitable harmonic oscillator bases jn^,, Uy), where n^, ny denote the numbers of oscillator quanta in the x- and 
y-directions. In the following we take the unit h = 1. The frequency ujq of the harmonic oscillator basis is determined 
so as to minimize trH in the adopted model space. The obtained values of loq are 7.51(fc = 0.2) and 7.13(fc — 0.6). 
The Hamiltonian matrix can be decomposed into submatrices due to the C4y symmetry. As in Ref. pO| we take up 
four classes of the one-dimensional representation which are labeled as ^j, ^2 , Bi, B2 according to their symmetry 
properties under reflection on the axes and diagonals in the x — y plane |lq| . (For instance, Ai is symmetric under 
both reflections.) The model space is spanned by the bases with < Ux + ny < 300, which gives the dimension of 
each submatrices as 5776, 5625, 5700, 5700. The diagonalization has been performed for each submatrices. Study of 
the nearest-neighbor spacing distribution confirms the character of the system suggested by the classical phase space 
structures, i.e., the Poisson like distribution for k = 0.2 and the Wigner distribution for k = 0.6 within each symmetry 
class. We also confirmed that the amplitude distribution of the wave functions for A;=0.6 show the Porter-Thomas 
distribution except for a singular peak at zero. 

In the following the results of the calculation will be shown for the states which belong to the symmetry class Ai . 
The results are similar for other symmetry classes. 

The basis state |7T.i,n2)ES belonging to the class Ai is written as 



\ni,n2)Es = v/ ^^^^^(|n-i,«2) + 1^2, ni)), (2) 

where ni and n2 are even integers and ni < n2. Because of the selection rule, the relevant part of the operator which 
contributes to the matrix element of the response should have the definite symmetry property. It should be noted 
that among the 5776 wave functions of the class Ai, those having very large energies are not quite reliable because 
of the limitation in the basis states. This is especially so when the response functions with large momentum transfer 
q are concerned. We may make an estimate for the range of validity by comparing the obtained level density with 
the semiclassical one. The comparison suggests that the maximum reliable energy to he E — 1000 ^ 1500 depending 
on the values of the parameter k in Eq. p). This maximum energy is contrasted with the largest energy eigenvalue 
E ~ 3000 obtained by the diagonalization. This limits the maximum value of the momentum transfer of the probe 
adopted below to be around q ~ 50, where the corresponding 'quasielastic peak' lies around ii^pcak = 9^/2 — 1250. 
This is confirmed by the calculation as shown later. 

III. RESPONSE FUNCTIONS 

We consider the response functions defined by 



w^^{n) ^ Y.\^jm)\H{n - {E, ~ E,)\ (3) 

where Q denotes a probing operator which connects the initial and the final eigenstates \i) and \j). In many cases of 
interest, the initial state \i) is set to the ground state of the system \gs) which belongs to the Ai symmetry class, in 
which case the index [i) is dropped. The response function shows the distribution of the state Q\i) over the energy 
eigenstates {|j)}. If, for instance, the initial state \i) has a simple structure as in the ground state of the harmonic 
oscillator, the state Q\i) and the response function will simply reflect the structure of the probe operator Q. We 
consider operators depending only on a single variable, say x, to see how the irregular behavior of the wave functions 
controlled by the parameter k may be reflected in the response function. One may rewrite the response function (H) 
in the form of the time-correlation function 

1 1"°° 
W^'^m - TT / ^^^'''* (Q(x(t))tQ(£(0))>„ (4) 
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where ( )i denotes the expectation value in the initial state \i). The probe Q is written as a function of the operator 

x{t) = e'"*x e"'-^* ^x+pj~ {x^ - kxf)f + • • • , (5) 

which is the solution of the Heisenberg equation of motion. In Eq. (0) we show also a short time expansion in terms 
of the operators at i = 0, e.g., a;(0) = x. A corresponding semiclassical expression for Eq. (Q) has been fully utilized 
in the analysis of Refs. llQ-flJ]. 

It is generally believed that the universal behavior of a dynamical system, i.e., if it is regular or chaotic, emerges 
in the fluctuation properties of the matrix elements of the operators, while their expectation values are strongly 
dependent on the specific dynamics of the system. Although the transition matrix elements for a chaotic system are 
known to generically follow the Porter-Thomas distribution, the energy-strength correlation such as the one contained 
in response functions is certainly dependent on the specific properties of the dynamics governed by the Hamiltonian. 
In this latter respect we note that the shape of the response function versus energy is constrained by a number of 
sum rules [g3[. Let us define 

/oo ^ 

dn n-w^^^n) = ^^(i?, -i?,)"l01QN)r (6) 
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for a given operator Q. The integer n may in general take negative values (usually for i — g.s. and j ^ g-s.), in 
which case the sum rule corresponds to the generalized susceptibility of the system. By increasing the value of n and 
subtracting the lower moments, e.g., in the form of the shifted moment or of the cumulant, one can regain finer and 
finer structure of the response function. Although one may recover the response functions by Mellin-transforming the 
sum rule values, the use of sum rules lies in the fact that in some cases the sum becomes a simple matrix element 
in the initial state. The latter may be calculated precisely and serves as a check for the accuracy of the calculation. 
The low n sum rules, in particular, sometimes become insensitive to the detailed dynamics and constrain the gross 
behavior of the response functions. 

As a probe Q of the response we first consider the operator i^ with an arbitrary initial state |z). We then fix 
\i) = \gs) and adopt the operator Qq = e"*^ which is closely related to an excitation of the system by an external 
probe characterized by momentum transfer q (and length scale 1/g). The operator x^ may be regarded as a long 
wavelength part of Qq, and is similar to the E2 operator of electromagnetic transitions. By changing the value of q 
in Qq, one can study in principle long as well as short distance structure of the matrix elements. 

In the actual calculation we consider only the operators symmetric under the reflection about x- and y-axes and 
also about the line x = y, i.e., 

Q = {Q)a, = \{Q{x) + Q(-x) + Qiv) -f Q(-y)), (7) 

where the arguments are explicitly written to show the dependence on the coordinates. For the initial state \i) in the 
class Ai this implies that in Eq. (^) only the states j belonging to the Ai symmetry class contribute. 

A. Response to the x^ Probe 

We first study the response function for the x^ probe. A typical example of the response to this probe is given in 
Fig.|l|. Here we show the response W'^^\Vl ^ E — Ei) for the z=500th initial state as a function of the energy E for 



k =0.2 and 0.6. For other initial states the main features are similar. We immediately see that the response consists 
of three clusters of strengths for both fe=0.2 and 0.6: Largest strength lies aX E = Ei with almost no other strength 
close to this peak, while two other clusters are located around E — Ei zL A, where A is slightly less than 2oJo, the 
expected value for a simple harmonic oscillator. For k — 0.2, the strengths are concentrated on a few states, while 
strengths are distributed over some energy range for k — 0.6. We now introduce creation and annihilation operators 
aj,, Qx, etc., of the oscillator quanta with frequency ujq, and decompose the operator x'^ as 

x^=Do + Dl + D2; Do=^{nx+ny + l), dI = ^{al^ + af). (8) 

One may then be tempted to assign the E — Ei peak to the response to the operator Dq, and two other clusters 
around E — Ei ± A to the operators D^ and D2, respectively. Figure (a) and (b) show the response to the operator 
Dq for k = 0.2 and 0.6, respectively. Although the state Do\i) is not proportional to \i), the fragmentation of the 
strength is restricted to only a few states in both cases. More important is the mixing in the state Dl\i). Figure g (c) 
and (d) show the response to the operator Dj for the i = 500th state. The distribution of strengths is considerably 
different between k — 0.2 and k ~ 0.6 cases. Strengths for k — 0.2 are seen to concentrate on a few states, while those 
for k = 0.6 are distributed over many states. 

The above features may be quantified by studying the number of principal components (NFC). The NFC for a 
normalized state \a) in terms of a complete set of orthonormalized states {\j)} is defined as: 

j 

The NFC becomes unity when the strengths are concentrated in a single eigenstate, while becomes A^tot when the 
strengths are equally distributed over the whole iVtot eigenstates. Figure ^ shows the NFC of the state Do\i) (with 
a suitable normalization) for each eigenstate \i) as a function of Ei, where the set {|j)} has been taken to be the 
eigenstate of the Hamiltonian. The NFC takes values 1 ~ 1.5 in most cases and decreases as the energy increases. 
General trend is not much different for k = 0.2 and k = 0.6. The situation is drastically different if we study the 
NFC for the state Dl\i) (again normalized) as seen in Fig. ^. For k — 0.2, the NFC remains small and does not show 
a marked energy dependence, while for k = 0.6, the NFC shows a rapid increase as a function of energy and takes a 
quite large value. We thus find that the D2{Dq) part of the operator x'^ is sensitive(insensitive) to the characteristic 
changes in the dynamics governed by the parameter k. 

There is another measure to see the difference between the two cases, k = 0.2 and k = 0.6, which can be obtained 
from the response function associated with x"^. In Fig. S the fraction of strengths (omitting the one for E — Ei) 
exhausted by two major states carrying largest strengths for each initial state \i) is plotted against Ei. For A:=0.2 
more than 60% of the total strengths is exhausted by the two major states and the distribution of the fraction is 
almost independent of the initial energy Ei, while for k = 0.6 they carry less than 50% and this fraction decreases 
as a function of the initial state energy. Since the dominant part of the strength associated with the operator Dq is 
contained in the initial state \i) and is omitted here. Fig. ^ shows the characteristics for the response to the operator 
D*2 (and D2) in accordance with the results from NFC. 

These studies imply that it depends strongly on the choice of the probe whether the difference in the character of the 
dynamics, namely regular or chaotic, may be easily seen in the response function. In the present case, the difference 
in the dynamics is not apparent for the probe Dq ^ x^ +p^, while it becomes quite significant for D2 ^ a:^ — p^ which 
is a probe, in a sense, 'orthogonal' to the unperturbed oscillator Hamiltonian ~ Dq. The response function for the 
probe x^ shows both characteristics. 

B. Response to the Probe Qq 

We now consider the response function for Qq = e"'^. We fix here the initial state to be the ground state. In this 
case the response is closely related to the situation of physical interest such as the inelastic electron scattering from 
the target in the ground state, where q gives the momentum transfer on the target. The symmetrized probe for Qq 
is given by Qq = (Q,)ai = ^e'"?* + e"''* + e'«« + e"'?*). The response functions W{q, n = E- Eg,) (for Qq) and 
W{q, n,) (for Qq) are calculated in terms of the elementary matrix element Qnmil) = {n\e^'^^\m) for a one-dimensional 
harmonic oscillator between the states with quanta n and m, which is given by 



Q„„(g) = ^"e^-z^^/-4^Ti"W (™>«), (10) 

y [n + a)] 



where a = to — n, z = q^ /2u>q and i"(z) denotes the associated Laguerre polynomiaL The functions Qnm{q) are 
calculated from recursion relations. 

We consider several q values corresponding to different resolution of the probe, the small q limit being related to 
the long-wavelength probe x'^ above. On the other hand, for large q values the operator resolves a fine structure of 
the system and the main strength of the response lies at high energies. If we use the short-time expansion in Eq. (1q) 
at this high 51 region, we can rewrite Eq. (m for the probe Qq using the Baker-Campbell-Hausdorff formula as 

where the dots denote operators with higher powers of t. The expression shows that the response is peaked at the 
quasielastic energy 9^/2 and has a width increasing with q and with, e.g., the momentum spread in the ground state. 
This holds precisely for a simple harmonic oscillator Hamiltonian, while in general is modified by anharmonicity 
effects. The limiting form of response at large q has been used to extract momentum distribution of complex system 
in terms of y— scaling analysis [E4|. In our case, as noted earlier, the model space of diagonalization limits the value of 
q around 50 with the corresponding limit ~ 1/q in the resolution of the wave function. This is much smaller than the 
length parameter l/y^ujQ of our oscillator basis. For the quartic oscillator the length scale will be modified from the 
simple oscillator value. One may define the characteristic length scale in the ground state by Xgs = {{gs\x^\gs)y^^. 
Calculated values of x~^ are 1.64 for k — 0.2 and 1.57 for k — 0.6. Thus the operator at, say q = 20, probes already 
a fairly fine structure of the system compared with the length scale of the ground state. The fact that the state 
Qq\gs) has an oscillation length scale 1/q also explains the occurrence of the quasielastic peak: The typical oscillation 
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length scale of the harmonic oscillator wave function at energy E ^ nujQ is ^J (xP-) jn ^ 1/ yE which becomes ^ 1/q 
in the region E ~ 9^/2. Thus, the state Qq\gs) will have the largest overlap with the states in the quasielastic region 
producing a peak in the response. 

Let us now consider the sum rules. Low n values of the sum Sn = 5„ {Qq) for the unsymmetrized probe are 
explicitly calculated to give Sq = 1, ^i = ig^, S2 — \q^{q'^ + f^-gs), etc., where the last sum rule is obtained from 

the virial theorem. For the symmetric probe the sum Sn = Sn {Qq) is not analytically obtained but is given by the 
expectation values as: 

5^1 = — q^(2 — cos2qx — cos2qy))„s, etc.. (12) 
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: for g->0, (13) 

: for q -^ 00, (14) 

where the values at g — *■ 00 are obtained under the assumption that the wavelength 1/q is much smaller than the 
typical length scale of the ground state wave function. These values are used to check the accuracy of the calculation 
within our model space, especially the one at large q which requires the matrix elements Eq.(nG) with large n. Note 
that these values are almost independent of fc, and the numerical calculation confirms that the fc-dependence is small 
indeed. It turned out that the limiting values (113) for the sum rule are satisfied already at q ~ 10. In view of Eq. 
(O) this result implies that the ground state expectation values of cosqx, etc., are almost zero, i.e., the resolution 
at (7 ~ 10 is already sufficiently fine for the ground state in accordance with the estimate given above. For higher n 
values the dependence of Sn on k is expected to become larger. Thus the gross structure of the response such as the 
total strength 5*0 and average energy Si/Sq is rather insensitive to the values of k. Detailed structure related to high 
n values of Sn should refiect the dynamics. 

Figure pi shows the response function W{q, ^) at q — 10, 30 and 50. The gross structures at a given q are similar 
for both k =0.2 and 0.6, and follow the behavior suggested earlier in this section: Not only the central energy follows 
f2 = q^/2 but the width of the response increases almost linearly with q. This should come out exactly from the sum 
rule if one had employed an unsymmetrized probe Qq. Fine structure of the response, on the other hand, is quite 
different for the two cases. For k = 0.2 the response has a rather simple regular structure: at q =30, for instance, the 
response is a superposition of a few structures with different sizes, each of which is centered around 9^/2 = 450 and 
is similar to the response of a harmonic oscillator given by 

Wi,^{q,n)^J2^(^~''''o)fn{^l fn{z)^^z"e-^. (15) 



In fact, by inspecting the wave functions one finds that these structures are related to the strong transition matrix 
elements of uncoupled (i.e., k = 0) quartic oscillators which are integrable. In contrast, for k — 0.6 the simple 
structure disappears and the values of the strength change drastically from one state to the other, although one can 
still see even at g = 50 a structure of the regularly spaced spikes as seen for k =0.2. 

For chaotic systems such as represented by random matrix theories, the amplitude distribution of wave functions 
or the matrix element distribution of an operator is known to follow the Porter-Thomas distribution. This may be 
contrasted to regular systems, where in many cases the quantum number imposes a selection rule of allowed transitions. 
Our result of response functions is in accordance with this generic behavior as far as the energy-strengths correlation 
is disregarded. In Fig. m we show the strength distribution of the response function for k =0.2 and 0.6 at q =30. For 
k =0.6 the distribution follows the Porter-Thomas form given by the dashed line and is quite different from that for 
fc=0.2. 

The question then arises: What is the nature of the persisting regular structure in the response functions of Fig.o 
for fc=0.6. This structure becomes more visible if one introduces a normalized response function defined by 

W„ormalizod(g,f^) = ~^. ' ^ p(»), (16) 

w{q,n) 

where W and p are respectively the response and the level density smoothed over energies. For the smoothing we 
employed the method of Strutinsky |g^ with the smoothing width of 20. This normalization procedure removes the 
gross structure effect of the response as constrained by the low order sum rules and enhances the embedded fine 
structure |22|. Figure || shows the normalized response for k =0.2 and 0.6. They show that the strengths in the 
regular spikes for k =0.2 are mostly redistributed for k =0.6 to produce smaller and smaller strengths to fill up the 
background, although one can still see the equidistant structure. The latter may be called an intermediate structure 
following Ref. |27| . This energy-strength correlation in the response function has been washed out in the strength 
distribution. 

The presence of an intermediate structure can be detected also in the response correlation function C(e) defined by 

C{e)= f dEW{q,E)W{q,E + €). (17) 



For equidistant structures such as the one for the free response (|15[), the correlation function gives again the regular 
pattern with the same spacing, e.g., 

^free(^) = E '^(^ " ("' " n)u;o) fn{z) fw (z) = J^ ^^^t^2.(i-cos^ot) ^ (18) 



with z = q^ /2ijjQ and /„ of Eq.(15). Figure g shows the correlation function for the response function at (7 = 30. In 
the actual calculation we used the normalized response function (Hq) in order to remove the gross shape effect and 
used the level number displacement 5i instead of the energy displacement e. The resultant correlation function was 
then smoothed with a smoothing width Ai = 4. We find the oscillator pattern arising from the intermediate structure 
in the response function for both k = 0.2 and 0.6. 

To understand why there still remains the intermediate structure even in the chaotic case, let us investigate the 
nature of the peak levels in some detail. We first pick up the peak levels carrying the largest strength among 
neighboring levels in the response function. For k — 0.2 the assignment of peak levels can be done without difficulty, 
while for k — 0.6 there may be an ambiguity. The qualitative results of the following analysis are however independent 
of this ambiguity. 

We first studied the NPC's of these peak levels in terms of the basis states (i.e., \j) = |ni,n2)ES in Eq.(||)), and 
found that they are markedly smaller than NPC's for other levels. This implies that the mixing of the basis states in 
these peak levels is smaller than other states. The nature of the peak levels may become clear from Fig. nfl, where 
we plot the quantity 

A^ = {i\\n^ - ny\/{n:, + ny)\i). (19) 

Here, black points correspond to the peak levels and bars indicate the average values over neighboring levels. For 
k = 0.2, the values of An for peak levels are nearly 0.9, almost twice of those for other levels. For k = 0.6, although 
it is not so evident as for k — 0.2, the values for peak levels are larger than the average. These facts strongly indicate 
that the peak levels are associated with the basis states of the type |0,to)eS: although considerably affected by the 
mixing with other states for k =0.6. In fact, these basis states are the only states excited by the probe Qq for a simple 
oscillator Hamiltonian. 



We may study the above results from the opposite direction. In Fig.|ll| we show the NPC of Eq.(^) for the basis 
state |a) = |ni,n2)ES by taking the eigenstates \i) of the Hamiltonian for the states \j) in Eq.(||). The abscissa is the 
basis number a, and the basis state type |0, m)ES is denoted by crosses. For both fc-values the basis states of the type 
|0,77i)es have the smallest NPC values, which shows that these basis states have the smallest spreading width caused 
by the mixing with other basis states. Thus the fact that there remains an intermediate structure in the response 
function may be restated in the doorway state picture ||2^: The probe Qq excites first the doorway states |0,r7i)ES 
which will then mix with other states causing the spreading of the strengths. As the spreading width for these specific 
states are smaller than the level spacing between the doorway states, the intermediate structure emerges even in the 
chaotic case. 

Since the states |0, m)ES may correspond to the classical isolated periodic orbits along the x and y axis, the relation 
between the intermediate structure in the response function and the scar M may be an interesting problem. Suppose 
that an initial wave packet \(f>{0)) is located at some point of the closed orbit having the period T. The wave packet 
will then semiclassically evolve along the closed orbit and will return to the initial position at each time interval T. 
Accordingly, the overlap of the wave packet at time t with the initial wave packet |((/>(0)|(/)(t))| will have peaks at 
t = nT, where n is an integer. The value of these peaks will decay due to the instability of the closed orbit like 
exp(— A/2t), where A denotes the Lyapunov exponent of the closed orbit M. Taking the state Qq\gs) as the initial 
wave packet 10(0)), the response function is nothing but the Fourier transform of the overlap |((^(0)|(/)(t))|. Therefore, 
the intermediate structure with peak level spacing D — 2t: /T and the spreading width of the peaks 7 = A may emerge 
in the response function if the condition 7/Z? < 1 is satisfied. In the present model, the period of the closed orbits 
along the x- and y- axes does not depend on the parameter k and is given by T = ^^^^(3)^ — ^.2A at _E = 1/2. 
Numerical calculation also shows that the Lyapunov exponents for these closed orbits sA E = 1/2 are A ~ 0.53, 
and 0.92 for k = 0.2 and 0.6, respectively. The values of 7/-D are then given by 0.44 and 0.77 for k = 0.2 and 0.6, 
respectively, both of them being smaller than unity. Note that j/D is independent of the energy. Thus, the existence 
of the intermediate structure may be explained also from the semiclassical point of view. 

IV. SUMMARY 

In this paper we studied properties of the response functions for a coupled quartic oscillator with several probes 
with a special attention to the difference between the regular and the chaotic cases. 

As a first example, we took the response to the probe x^. Since the response function is determined by the probe 
as well as the nature of the wave functions of the system, we must pay attention also to the character of the probe. 
For instance, the operators for which the diagonal matrix element becomes a main coniponent are not adequate to 
see the difference of the dynamics. For the operator x'^, it can be decomposed like Eq. (jg), and the operator Do has 
such a character, while for the operators D^ and D2 the non-diagonal matrix elements are important. Therefore, by 
removing the strengths associated with the operator Dq we can see the difference of the response between the chaotic 
and the regular cases, namely more spreading of strengths for the chaotic case, which can be quantified with the 
number of principal components. 

Next, we considered the response to the probe Qq — e"^^. The response function at a given momentum transfer q 
is related to the time-correlation function of the operator with a resolution l/q in the coordinate space. It was shown 
that the gross structure of the response function is similar for the chaotic and the regular cases as constrained by 
global sum rules. On the other hand, the difference is reflected on the fluctuation, as seen in the strength distribution 
( the histogram of strengths) . Moreover, we detected the intermediate structure (i.e., typical energy scales) even in the 
chaotic case which can not be expected for the random-matrix model. We found that the existence of the intermediate 
structure is due to the fact that the spreading width of the doorway states is smaller than the level spacing of the 
doorway states, and also indicated its relation to the scar phenomenon. It would be interesting to study this structure 
from a different point of view, e.g., the semiclassical theory of responses ]l0[-fi3[ based on the periodic orbits. 

In this paper for the sake of simplicity we restricted the discussion to the transitions between the states belonging to 
the same symmetry class. In the realistic situation, however, the transition connecting states with different symmetry 
classes may also occur at same time. It is well known that the level spacing statistics drastically changes when we 
consider the levels belonging to different symmetry classes simultaneously. Thus, it is also interesting to see what 
happens for the response function to the probe connecting different symmetry classes. 

The authors thank M.Matuso for valuable discussions. They thank also P.Schuck for a discussion about the 
semiclassical description of response functions. 
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FIG. 1. Response to the probe x^ for the initial state \i = 500) as a function of the energy E for k =0.2 (a) and 0.6 (b). 
Contained levels in the displayed energy regions are 459th to 542th for k = 0.2 and 462th to 554th for k — 0.6. 

FIG. 2. Response to the operator Do for the initial state ji — 500) as a function of energy E for k =0.2 (a) and 0.6 (b), and 
the one for the operator D^ for k = 0.2 (c) and 0.6 (d). 

FIG. 3. The NPC NpJ for the state Do\i) as a function of the energy of the initial eigenstate Ei for k = 0.2 (a) and 0.6 (b). 

FIG. 4. Same as Fig. | but for D^. 

FIG. 5. The fraction of the x^ strengths (omitting the one for E = Ei) carried by the two major states to the total strengths 
is plotted for each initial state \i) as a function of the initial state energy Ei. 



FIG. 6. Response function W{q, E) for ifc = 0.2 at g = 10 (a), 30 (c), and 50 (e) and for fc = 0.6 at g = 10 (b), 30 (d), and 
50 (f). Note the changes in the scale of the vertical axis. 

FIG. 7. Strength distribution P{S'^^'^) of the matrix element 5" = \{j\Qq\gs)\^ at g = 30 for k = 0.2 (a) and 0.6 (b). Strengths 
are normalized as Eq. (lid). Dashed line shows the Porter- Thomas distribution. 

FIG. 8. The normalized response function Eq. (Ilfl) at q =30 for k =0.2 (a) and 0.6 (b). Compare with the responses at 
q =30 shown in Fig. |6[ 

FIG. 9. Smoothed correlation function C{5i) for the normalized response function to the probe Qq with g = 30 as a function 
of the level number displacement Si for k = 0.2 (a) and 0.6 (b). Smoothing width Ai = 4 is adopted. 

FIG. 10. The quantity A„ for k = 0.2 (a) and 0.6 (b). Black points denote the peak levels for the probe Qq with g = 30 
and bars show the average over neighboring levels. 

FIG. 11. The NPC A^pc for basis state \a} belonging to the symmetry class Ai for k = 0.2 (a) and 0.6 (b). Horizontal axis 
shows the level number of the basis state |q). Cross points correspond to the basis states type |0, m)ES. 
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